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SyYNopPsIS 


Analytical expressions for the effective width of cylindrical shells are pre- 
sented herein. These formulas are used for the analysis of a cylindrical shell 
roof acted upon by horizontal longitudinal forces. Also studied are the influ- 
ences of foundation movements, temperature changes, and shrinkage. 

A comparison is made between theoretical results and experimental meas- 
urements obtained from a model of an actual cylindrical shell roof constructed 
in the approximate ratio of | to 30. It is established that there is satisfactory 
agreement between the theoretical analysis-and the test results. 


INTRODUCTION 


It is a common misconception that the use of shells as structural elements 
is a development of the twentieth century. Two thousand years before, the 
Romans had constructed the Pantheon, which was covered by a spherical shell 
with a diameter of approximately 145 ft. The cupolas of the cathedral in 
Florence (Italy) and St. Peter’s Basilica in Rome (Italy), built during the 
fifteenth and sixteenth centuries, respectively, and having a diameter of ap- 
proximately 140 ft, are examples of remarkable craftmanship, based only on 
experience and intuition. The contribution of the twentieth century to this 
type of construction consists in the development of a rational analysis which 
can be used to predict the state of stress in the shell. 

The development of the shell theory was initiated by A. E. H. Love (1)* 


+? Numerals in parentheses, thus (1), refer to corresponding items in the Bibliography (see Appendix 


Il). 

who derived the differential equations for curved plates. E. Meissner (2) (3) 
succeeded in integrating these equations for spherical and conical shells. F. 
Bauersfeld (4), U. Finsterwalder (5), and F. Dischinger (6) made possible the 
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application of the theory to the analysis of actual structures. A number of 
reinforced concrete shells were built in Europe immediately following World 
War I, and there were concurrent theoretical contributions. 

The arena in Hershey, Pa., built in 1936 (7), was the first modern cylindrical 
shell roof constructed in the United States. In the years that followed, and 
especially during and subsequent to World War II, a number of shell roofs 
were built—to be used as storage houses, factory buildings, airplane hangars, 
sport arenas, and armories—covering a total of 10,000,000 sq ft. 

At present (1953) the largest span for a cylindrical shell roof is 340 ft (the 
hangars at Rapid City, S. Dak. (8) and at Limestone, Me. (9)). These large- 
span structures accentuated certain problems which were of secondary impor- 
tance in the previously built smaller structures. Because there is a tendency 
toward increasing the span of shell structures, a careful analysis of the pertinent 
problems is mandatory. 

The supporting ribs in cylindrical shell roofs can be designed as T-sections 
to include the additional effective section provided by the shell. A study 
(10) of the interaction between the ribs and shell shows that an “effective 
width” of the shell can be determined and the combined section used in a rela- 
tively simple arch analysis of the rib structure. 

Lateral horizontal forces, caused by wind action on the front and back 
doors of large-span cylindrical shell roofs, can reach a considerable magnitude. 
The influence of these forces on the strueture cannot be considered secondary. 
A difference of temperature inside and outside the building can produce rela- 
tively high thermal stresses. Shrinkage and plastic flow greatly affect the 
stability of the structure. These are a few of the problems to be considered 
in the analysis of large shell roofs. 

During the period from January 1, 1949, to February 15, 1951, a research 
program on shell arch roofs was conducted at the Fritz Engineering Laboratory 
of Lehigh University at Bethlehem, Pa., with special emphasis on the pre- 
viously cited problems. A careful theoretica! study of the problems was made, 
and the results were checked experimentally on the model of a cylindrical shell 
roof shown in Fig. 1. This paper presents three contributions to the analysis 
of shell roofs made during the course of the research program. It will be 
shown that the model tests were in excellent agreement with the theory. 


ErrectivE WiptH oF CYLINDRICAL SHELLS 


The stress distribution in the flange of a T-beam with a straight axis, sub- 
jected to bending in the plane of the rib, is not constant over the width ef the 
flange. This condition is shown in Fig. 2. Simple beam theory, based on the 
Navier-Bernoulli hypothesis (a plane cross section remains plane after bending) 
can be utilized if the actual width of the flange is replaced by the effective 
width b. T. von Karman, Hon. M. ASCE (11), was the first investigator to 
derive the correct theoretical expression for this effective width (12). Shear 
lag is responsible for the considerable damping of the direct stresses in the 
flange in the direction of the rib. 

In Fig. 3 there is shown a cylindrical shell stiffened by a rib in the circum- 
ferential direction. Arbitrary loads, acting in the plane of the rib, will cause 
direct forces Ng per unit width in the circumferential direction. The actual 
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cylinder can be replaced by a ring of width b for which a constant direct force 
(N¢)z—0 (equal to (N4),~. for the cylinder) over the entire width is assumed. 
By determining the width b so that 


the stresses in the ribs are obviously identical for both structures under an 
equal load system. If width b is known, determination of the rib stresses 
becomes a problem that can be solved by the use of the simple beam theory. 

Analytical expressions for the effective width, based on the general bending 
theory of cylindrical shells (13), have been derived by one of the writers (14). 
There is also a short report on this study in which the analytical values for 6 
are compared with experimental results (10). 

C. B. Biezeno and J. J. Koch (15) have derived expressions for the effective 
width of a shell that extends indefinitely on either side of a rib. This is a 
special case (Bl = © in Fig. 4) of the general solution in which the rib can be 
any distance / from a free boundary. It is shown that, essentially, two effects 
goverrt the effective width of cylindrical shells: (a) The lag of the direct shear 
forces in the plane of the shell (as in the case of a T-heam with a straight axis) 
and (6) the radial escaping of the shell under the circumferential direct forces 
Ng. In general, 6 can be written (14) as 


in which r is the radius of the eylinder, t denotes the thickness of the cylinder, 
and K equals f(A, Bl) (Fig. 4). 

The parameter » (Fig. 4) takes into account the influence of the direct 
shear forces on the effective width 


in which n depends on the variation of the stresses in circumferential direction 
and is the number of complete waves made by a harmonic function around 
the circumference of the cylinder (Fig. 4). For equal to zero, the stress 
distribution has axial symmetry, and the effective width has its maximum 
value. The parameter 8 / is dependent on the length of the overhang of the 
shell / and on the shell constant (13a) 


As shown in Fig. 4, the shell is assumed to. extend infinitely on the left side. 
This can be assumed safely if the product Bs exceeds 2.4, in which s is the 
distance of the left end from the rib (14a). 

Taking the effective width 6 as the flange and the rib as the web of a T- 
section, the Bernoulli-Navier hypothesis leads to the correct values for the rib 
stresses. Also, denoting the circumferential stress along the connecting line 
of rib and shell by oa, the total direct force F in the shell.is the integral of all 
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direct forces (Ng = og t) over the entire length of the shell or is equal to the 
product of the stress oa, the thickness ¢, and the effective width b. Thus, 


For a given value of F and a given set of boundary conditions, all forces and 
moments in the shell can be determined. Tables and graphs which reduce 
these computations to a minimum are available (14). 

Knowledge of the effective width of « cylindrical shell simplifies an involved 
problem of elasticity to a commen problem of the ordinary beam theovy. . This 
effective width will be used subsequently in the analysis of shell roofs. 


Tue Errecr oF Horizonrat Forces 

The action of the wind on the front and back doors of cylindrical shell roofs 
produces considerable horizontal forces. This type of loading will be analyzed, 
and numerical results will then be compared to test results obtained on a model. 

During construction, the roof is erected in units in order to re-use the forms 
and to provide expansion joints. One unit of a cylindrical shell roof is shown 
in Fig. 5 and it is subjected to horizontal forces in the direction of the axis. 
On a cross section a shear force V, a bending moment Wx, and a twisting 
moment M 7 are acting. It is necessary to find the stress distribution resulting 
from these forces and moments. 

In Fig. 6 there is shown a part of a cylindrical shell of infinitesimal width 
rdg. The forces and moments acting on the shell element are defined in Fig. 7. 
The direct force N, in the axial direction is equal to the applied uniformly 
distributed load p at the rib xz = 21. The force V, equals zero at the rib x = 0. 
Assuming a linear variation of N, along the z-axis, then 


| 
The assumption of a linear variation of N, is especially true for shell roofs whose 
width is relatively small in comparison to the span L. The equilibrium of an 
infinitesimal shell element dz r d@ in the z-direction results in (Fig. 6) * | 
ON, ON eos 


Replacing 0N,/0z by its value from Eq. 66 results in 


ON 
= ‘ 


By integrating Eq. 7b and considering the fact that the shear force Nig is zero 
at the center of the span, because of symmetry, \4: becomes 
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It follows that Ny, is constant for a section in which @¢ is a constant and 
varies in proportion to the angle @ along the span of the shell. From the 
equality of the shearing stresses (rg: = 729) the direct shear forces, Nzg, in a 
radial direction can be considered equal to the direct shear forces Ng:. Thus, 


Nig = Nez 


Eq. 9 is a reasonable approximation, but it should be realized that Nig is not 
exactly equal to Ng; as a result of the curvature of the shell (16). 

If imaginary support forces Nz, acting along the outer ribs (x = 0 and 
az = 21) are assumed, the shell is in equililrium under the load p and is sub- 
jected only to shear forces Ng and Ng, and to direct forces N,. Actually 
these imaginary forces are not acting. Introducing the shear forces —Nig 
along the connecting line of the shell and the ribs, the imaginary forces Nig 
are eliminated. This procedure is similar to the method of moment distribu- 
tion. In that method, the joints are locked by imaginary moments which 
afterward are eliminated by introducing the moments in the opposite direction. 
The interaction between the rib and the shell is accounted for by considering 
the effective width of the shell as a flange of the rib. This effective section, 
consisting of the rib and the effective width, is used in analyzing the arches. 

Analysis of the Ribs.—The rib at x = 21 is acted upon by shear loads, Nz¢, 
per unit length, acting along the connecting line of the rib and the shell. The 
loads are considered positive as shown in Fig. 8. On the rib z equal to zero, 
the positive N,4-forces are acting in the opposite sense. Fig. 6 shows the sign 
convention for the forces Ng. The distance from the loads Ng to the cen- 
troidal axis of the effective section is equal to zg. The rib is an elastically re- 
strained arch, having only two redundants as a result of the symmetry of the 
structure and the load with respect to the center line. In Fig. 8 one half of one 
rib is shown with the horizontal thrust H, and the bending moment M, as 
redundants at the center. The two geometrical conditions furnishing the two 
equations for the determination of H, and M, are a known horizontal move- 
ment 6x of the foundation and a known rotation 0x. For a fully restrained rib 
the conditions are 


In the case of an elastically restrained arch, with known measured founda- 
tion movements, the conditions are 


ox (5k) mens 
and 
Ox = (Ox) meas Mx 


The terms (8x) meas 2Nd (OK) mens Aare the measured horizontal displacement 
and rotation of the abutment, respectively. The symbol x is the coefficient of 
elastic restraint, and the product —«x Mx, in which Mx is the restraining mo- 
ment, represents the elastic end rotation of the rib. 
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In general, the shape of the center line of the rib and the variation of the 
moment of inertia J of the effective section require a numerical procedure for 
the determination of the two redundants H, and M,. For the special case of 
a circular center line and a constant moment of inertia, the computation for 
H,and M, can be accomplished by direct integration. Extension can be made 
to other cases such as a parabolic shape of the rib and J equal to J,/cos ¢, or 
to a numerical solution. 

For the statically determinate system shown in Fig..8, the normal force NV, 
and the bending moment M,, caused by the shear load Nz4, can be obtained 
by integrating the contribution of the distributed shear load. If @ is the angle 
at which N, and M, are determined, then N,4 rdw represents the shear load 
applied over an infinitesimal rib length r dw at any angle w between zero and ¢. 
Then from Fig. 8 (0), 


N,= -f cos (@ — w) dw.............. (12a) 


Substituting from Eq. 9 and integrating Eq. 12a results in 
2 
N.= 
Similarly, the contribution to the moment M, of the incremental shear load 


is Nig r dw} Z, — r-[1 — cos (@ — w)]}, in which r, = r + 2q and is the radius 
of the effective section (Fig. 8(b)). Therefore, 


M,= Nig riza — r-[1 — cos (@ — w)]} dw........ (13a) 
Substituting from Eq. 9 and integrating Eq. 13a results in 


The total normal force N and the total bending moment M are (Fig. 8(a)) 


and 
M = M, + H.r. (1 — cos¢) + M,.............. (45) 
In order to compute H, and M,, the horizontai displacement 6x and the rota- 
tion of the abutment 6x must be determined. Choosing a virtual load system 
as shown in Fig. 8(c), P’ equal to unity at the abutment produces a virtual 
normal force N’ and a virtual bending moment M’ at an angle ¢ in the rib. 
These qualities are related by 


3 
> 
3 
' 
(16) 
M’ = —r, (cos @ — cos 
| 
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By use of the work equation, the displacement 6« of the actual system is ex- 
pressed as the work done by the virtual load (P’ = 1): 


as 
El 

in which 
ds = r,d@ 


and A is the cross-sectional area. The influence of the normal force N on the 
displacement is extremely small and therefore can be safely disregarded. 
Assaming that the product EF J is a constant and replacing M’ and M by their 
values from Eqs. 17 and 15, respectively, Eq. 18a becomes 


1 
6x = aif (cos — cos Ox) 
x [M, + H.r,. (1 — cos + M. Jr. dg. . (19) 
Replacing M, by its value from Eq. 136 and integrating Eq. 19, 


— in = Me (sin ox — ox bx) + Her, (sin ox + sin 2 ox 


e 


2 


3 dx ox cos bx) + {(: cos $x 


2 
sin dx — 24x + + 008 


The rotation 60x of the abutment is determined by the same procedure. The 
virtual moment M’x equal to unity of Fig. 8(c) produces a bending moment 
M’ equal to unity for any angle @ such that 


M’'k 


Inserting Eqs. 21 and 15 into the work equation, 


(Mas 


results in 
6x = El H.r. (1 — cos¢) + 


If M, is substituted by its value from Eq. 136, the integration of Eq. 23 leads to 


= + Here (ox — sin x) 


2 3. 
[( — ox + sin 
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Eqs. 20 and 24 each have the two unknowns, H, and M,.. The normal force N { 
and the bending moment M for any angle ¢ are determined from Eqs. 14 and 
15 following the solution for the redundants. The stresses in the ribs are 
determined from 


N Me 


in which A and / are, respectively, the cross-sectional area and the moment of 
inertia of the effective section, consisting of the rib as a web and the effective 
width of the shell as a flange. The symbol z represents the distance of a fiber 
from the centroical axis of the effective section. 

Coefficient of Elastic Restraint —In Eq. 11 a coefficient of elastic restraint x 
was introduced. In arch bridges this coefficient provides for the elasticity of 
the soil and the abutment. In general, the arch cannot be considered fully 
restrained at the springing line. j 
Dtring the tests performed on the model shown in Fig. 1 another consid- 


eration led to the introduction of x. The section A~A in Fig. 1 gives a typical 
construction detail at the springing line of the rib and the shell. The latter is 
supported by a relatively flexible edge member. The rib, however, terminates 
in a very heavy end wall and can be regarded to be fully restrained. There q 
are two extreme cases that can be considered: (1) The shell is held rigidly 
along the edge member and (2) the shell has a free edge. In the first case the 
N4-forces of the shell will be resisted by the support, and the effective width 
is sensibly constant down to the edge member as in Fig. 9(a). If there is no j 
support at all, the effective width must be reduced to zero at the springing line 
as in Fig. 9(b). The actual condition lies between these two cases. Since : 
this disturbance of the effective width of the shell has only a local effect, the | 
reduction of the moment of inertia of the effective section can be accounted 
for by considering the value of J as being constant to the springing line and 
by assuming an elastic restraint for the rib. The higher flexibility of the 
effective section in the end zone is therefore concentrated at the springing line. | 


It is obvious that the magnitude of the coefficient of elastic restraint depends 
on the effectiveness of the edge member in supporting the shell. No theoretical 
analysis was made of x during the investigations. 

Determination of the Shell Forces—The stress along the connecting line of f 
the rib and shell o, is computed by use of Eq. 25, following the determination | 
of the normal force N and the bending moment M of the effective section. 


The total circumferential direct force F of the shell is given by Eq. 5. 

The problem consists of finding the forces and moments in the-shell for a ‘ 
given F and a given set of boundary conditions. In the present case the | 
boundary conditions for the shell at the ribs x = 0 and x = 2/ are 


in which M, is the bending moment of the shell in the axial direction. The 
condition given in Eq. 26 presupposes the neglect of the torsional stiffness of { 
the rib. The influence of the middle rib on the stress distribution can be | 
disregarded, the forces and moments being rapidly damped with increasing 

distance from the end ribs. The actual computations for these forces and | 
moments are readily made by use of tables (14). 
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Deflections of the Rib.—Knowing the area A and the moment of inertia J of 
the effective section, the normal force V, and the bending moment M for any 
angle @ along the rib, the deflection of the rib is readily determined from the 
work equation. Placing an imaginary load P’ equal to unity at the point at 
which the deflection is desired, the work of this unit load as a result of the 
deformations caused by the actual load system is equal to the deflection under 


consideration, or ° 
M’ M ds N’Nds 


in which N’ and M’ are caused by the imaginary load, and M and N are a 
result of the actual loads. 

Numerical Example.—The foregoing theory will be applied to the analytical 
solution of a model of a shell roof. Fig. 1 gives the dimensions of the model 
of the hangar at Rapid City built to the scale of 1 to 30. Structural steel was 
chosen as material for the model, the dimensions (a shell thickness of approxi- 
mately } in.) leaving no other choice. The dimensions and properties not 
shown in Fig. 1 are as follows: 


Item Value 
30 X 10° lb per sq in. 
Coefficient of elastic restraint, 0.2 
The measured movements of the supports were 
5x, in inches $x, in radians 
Rib x 102 x 103 


Effective Section of the Outer Ribs.—The effective width of the shell is first 
computed. Fig. 4 can be used to determine the effective width under general 
conditions. The coefficient 8 is found, from Eq. 4, to be 


vrt 


To determine the number n, the stress distribution along the rib in a cir- 
cumferential direction must be known. Anticipating the final results, Fig. 10 
shows the stresses in the ribs. The stress o, (oq has the same variation) in the 
lower fiber has approximately the form of a cosine function, with a half-wave 
length of from 0.35 radian to 0.40 radian. Thus, 


= 
= half-wave length — 2 ox (29) 


The right side of Eq. 29 is an approximate expression for n in terms of the 
angle @x (shown in Fig. 8) and holds for this particular case only. Another 
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load would produce another variation of the stress oy. Substituting the 
numerical value for @x into Eq. 29 results in n equal to 8.03. The coefficient 
A is found, from Eq. 3, to equal 0.26. 

The two outer ribs have no overhang, and the distance to the middle rib 
can be considered to be infinite. Entering the chart in Fig. 4 with the values 
of B l equal to zero and X equal to 0.2¢, the coefficient K is found to equal 0.38. 
It can be seen that for values of X less than 6.75 the coefficient K and the 
effective width are sensibly constant. The effective width b, from Eq. 2, 
becomes equal to 1.3€ in. 

The area A and the moment of inertia J of the effective section are then 
computed. In Fig. 11 (a) there is shown the pertinent information for the outer 
rib. Further data for the outer rib are listed: 


Value 


Moment of inertia, J............... 0.536 in.* 


The ratio of the moment of inertia of the effective section to the moment 
of inertia of the rib is equal to 1.35. Thus, the shell increases the bending 
stiffness of the rib by 35%. The shear loads Nz4 acting on the ribs are deter- 
mined from Eq. 9. 

All values are available to solve Eqs. 20 and 24 for the two redundants 
H,.and M,. The solution is made for the assumption of a rigid restraint and 
an elastic restraint with measured displacement and rotation of the foundation. 
The following values are determined: 


Restraint H-, in pounds Mc, in inch-pounds 
— 1,457 +4,371 


On the basis of these values of the redundants, the direct stresses in the 
ribs were determined and are shown in Fig. 10. 

The computation of the direct forces Ng and the cross bending moment M, 
in the shell followed the procedure previously outlined. The results are pre- 
sented in Fig. 12. 

Experimental Investigation and Comparison—An extensive experimental 
study on the model shown in Figs. 1 and 13 was made. Strains were recorded 
by 44 rosette (type AR-1), 137 cross (type AX-5), and 81 single (type A-5) SR-4 
electrical strain gages. Displacements were measured by 20 Ames dial gages 
(with an accuracy of 0.001 in.), and the rotations were checked by two level 
bars (with an accuracy of 0.0002 radian). The actual loading consisted of ten 
equally spaced horizontal single loads as shown in Fig. 5(b). By virtue of 
St. Venant’s principle the differences between the theoretical loading (assumed 
to be uniformly distributed) and the actual loading can cause only local differ- 
ences in behavior. 

The stresses in the ribs computed from the strain-gage readings are shown 
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Fig. 12.—Forces anp MoMENTS IN THE SHELL 
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in Fig. 10. There is good agreement between the test results and the theoreti- 
‘al analysis if the analysis is made with an allowance for the elastic restraint 
of the rib. The introduction of «x is fully justified. In Fig. 12 the experi- F 
mental Ny (the direct force per unit width of the shell in the circumferential 
direction) and M, (the bending moment per unit width in the axial direction) 
of the shell are plotted. 

The agreement between the analytical and the theoretical results is quite 
close. It is of interest to note that the stress distribution is perfectly anti- 
symmetric with respect to the middle rib. The middle rib is unstressed and 
therefore does not add to the strength of the structure under the present q 
loading condition. The rib at x equal to 21 has exactly the opposite stresses q 
of the rib at x equal to zero. This fact is also demonstrated in Fig. 14 in which 
the measured vertical deflections of both outer ribs are plotted. 


Conclusions Concerning Horizontal Longitudinal Forces.— 


1. The analysis of the effect of horizontal lateral wind loads, on cylindrical 
shell roofs was fully confirmed by the experimental investigations on the model. 

2. For actual shell roofs of the type constructed at Rapid City, the wind 
pressure on the front door (as prescribed by building codes) is from 30 lb per 
sq ft to 35 lb per sq ft, or approximately 1000 lb per lin ft of the rib. This q 
results in maximum rib stresses of approximately 400 lb per sq in. Obviously, 
a stress of this magnitude cannot be considered a secondary one. By com- 
bining a number of shell units, these stresses can be lowered considerably. 


FouNDATION MOVEMENTS 


The horizontal thrust of long-span shell roofs are of sufficient magnitude ' 
to offer serious foundation problems. In certain cases, tension ties have been 
found necessary to balance the horizontal thrust acting on the abutments (8a). 

It is usually impossible to prevent any foundation movements, except by , 


special devices such as artificially stressing the tension ties between the abut- \ 
ments. 

The general procedure for the computation of the forces in arches subjected | 
to foundation movements can be used if the interaction, occurring between the i] 


ribs and the shell, is taken into account. This will be done by taking as a 
cross section of the rib the effective cross section formed by the rib and the 


effective width of the shell. ' 5 
Analysis of the Ribs —The horizontal displacement 6x and the symmetrical i 
rotations 0x of the two abutments are considered. Unequal rotations @«; and a 
6x, of the left and the right abutments can be solved by superimposing sym- a 
metrical rotations, i} 

4 

{ 

and antisymmetrical rotations, 

The advantage of this procedure is to reduce the number of redundants from ‘ 
three to two. {| 
The rib, shown in Fig. 15(a), has two redundants—the horizontal thrust . 
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Hx and the end moment Mx. Introducing the virtual load H’x equal to 
unity and the virtual moment M’x equal to unity, the horizontal displacement 
5x and the end rotation 6x, respectively, can be determined by use of the work 
equation. 

The actual normal force N and bending moment M (Fig. 15(a)) are 


and 
M = Mk — Hxr,. (cos @ — coS Ox)... (33) 


Because H’x is equal to unity (Fig. 15(c)), 


N’ = 
and 
M’ = —r, (cos @ — cos ox) 


M’ M ds N'’N ad 


eK 


Hence, 


Assuming the products E / and E A to be constant, inserting Eqs. 32, 33, 
34, and 35 into Eq. 36, and integrating leads to 


6x = Mx (ox cos — sin gx) + G+ cos” 
I 


— 2sin2 + Ar (3 dx + sin2 ox)| (37) 
€ 

If M’x and M’ are made equal to unity and N’ is made equal to zero, the 
work equation provides the end rotation 6x. The reduction of the effective 
width in the edge-member zone is accounted for by a coefficient of elastic 
restraint x. Thus the product of the external moment M’x equal to unity of 
the virtual load system times the rotation (—x Mx) caused by the elastic 
restraint under the actual load system gives a contribution to the external 


work: 
M’ M d: N di 


Substitution of Eqs. 32 and 33 into Eq. 38 and integrating yields 


in which 


For any given values of 6x and @x, the two redundants Hx and Mx are deter- 
mined by Eqs. 37 and 39. The direct stresses in the ribs are determined by 
use of Eq. 25. 


| 
i 
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Value of @ by 


Fic. AND STRESSES IN THE MODEL As A 
or FounDAaTION MoveMENTS 


Determination of the Shell Forces.—The procedure for determining the shell 
forces resulting from foundation movements follows the steps used in the case 
of horizontal longitudinal forces. For the two outer ribs and the middle rib 
the corresponding F-forces are computed from Eq. 5 following the determina- 
tion of the stress oq along the connecting line of the rib and shell (from Eq. 25). 
Knowing the boundary conditions for the shell and the F-forces, the stress 
distribution in the shell can be determined. This analysis can be found else- 
where (14). 

Numerical Example-——The model shown in Fig. 1 is analyzed for a given 
horizontal displacement 6x and end rotation 0x. The effective width of the 
ribs must first be computed. In order to know the number n, used in Fig. 4 
to determine the effective width, the stress distribution along the rib must be 
known. Fig. 16 shows that, as a first approximation, the variation of the 
stress a, in the lower fiber of the middle and outer ribs is a cosine function 
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with a half-wave length of approximately 0.75 radian. Hence, 


This value of n is a first approximation if symmetrical foundation movements 
are considered. Inserting the value of ¢x (0.5866 radian) into Eq. 41 results 
in n being equal to 4.14. From Eq. 3, A is found to equal 0.14. Thus, for 
the middle rib 6 / is equal to 4.42 (greater than 2.4) and for the outer ribs Bl 
is equal to zero. 

From Fig. 4 it is obvious that the values of K for (a) \ = 0, Bl = © in 
case of the middle rib, and (b) X = 0, Bl = 0 in case of the outer ribs are 
sufficiently accurate for the determination of the effective width. Thus, for 
the middle rib K is equal to 1.52 in. and b is equal to 5.43. in. For the outer 
rib K is equal to 0.38 in. and 6 equals 1.36 in. 

The following is a list of values for the effective section of the middle rib. 


Item Value 
1.961 sq in. 
— 1.382 in. 
0.731 in. 
Distance, Za 
Radius, 108.731 in. 
Moment of inertia, 7 0.921 in.‘ 


The ratio of the moment of inertia of the effective section to the moment 
of inertia of the rib is equal to 1.84. The considerable increase (84%) of the 
bending stiffness of the middle rib as a result of including the effective width 
of the shell is especially noteworthy. 

On the model a horizontal displacement, 6x, of 0.3300 in. and an end rota- 


tion, dx, of 7.753 K 107‘ radian were induced. 
Assuming for x the value, 


and inserting into Eqs. 37 and 39 the proper values, the horizontal thrust Hx 
and the end moment Mx are determined as follows: 


Rib Hk, in pounds Mx, in inch-pounds 
— 17,850 
— 10,350 


It should be noted that x depends on the type of loading; for the case of 
horizontal lateral loads a different value for x was used. 

The normal force N and the bending moment M for an arbitrary angle @ 
along the rib are given by Eqs. 32 and 33. The rib stresses (determined from 
Eq. 25) are plotted in Fig. 16 for the middle and outer ribs. 

Fig. 17 shows the normal force in a circumferential direction Ny, and the 
bending moment in axial direction M, of the shell. 

Experimental Investigation and Comparison.—The same model (Fig. 1), as 
used for the test of lateral horizontal loads, was subjected to a horizontal 
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The measured displacement 


foundation displacement as shown in Fig. 18. 
and the end rotations are as previously given. 

The stresses in the ribs and the NVy-force and M,-moment computed on the 
basis of the SR-4 strain-gage readings are plotted in Figs. 16 and 17, respec- 
tively. The theoretical analysis and the test results again compare favorably. 
The discrepancy of a few experimental M,-values in Fig. 17 does not influence 
the over-all agreement. 


Conclusions Concerning Foundation Movements.— 


1. By taking as the effective width of the shell the flange of a T-section 
(the web of which is the rib), the analysis of the ribs for foundation movements 
“an be accomplished by use of the simple arch theory. Following the deter- 
mination of the stresses in the ribs, the stress distribution in the shell can also 
be computed. Experimental results on the model support the theoretical 
analysis. 

2. Eqs. 37 and 39 show that, for given foundation movements 6x and @¢x, 
the moment Mx and the horizontal thrust Hx are proportional to the moment 
of inertia J of the effective section of the rib. Neglecting the influence of the 
normal force on the deformations, which is completely negligible, 


(43a) 


and 


in which fi(@) and fo(@) are functions of ¢. 
The stress in the rib becomes, from Eq. 25, 


A I A 


Thus, ¢@ is approximately proportional to the distance z of the fiber from the 
centroidal axis of the cross section because the first term of the right side of 
Eq. 44 is relatively small. The importance of these derivations will be ex- 
plained subsequently. 


Some SpectaL 
Uniform Temperature Change and Shrinkage of the Concrete-—-Suppose a 
shell roof undergoes a uniform temperature change of At®. If the whole strue- 
ture would be supported as a simple beam, this temperature change would 
sause an increase AL of the span L: 


in which @ is the coefficient of thermal expansion. Actually this increase (AL) 
is impossible because the abutments are restrained. Therefore the original 
span L must be restored by diminishing the new span (L + AL) by AL. The 
latter is equal to a horizontal displacement, 6x, previously investigated under 
the heading, “Foundation Movements.” By substituting into Eqs. 37 and 39 


the values, 


and 
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Fig. 18.—Lasoratory Movev Beine Sussectrep To HorizontaL FounpDATION DisPLACEMENT 
0.625 In. 0.118 in. 0.625 In. 


0.625 In. 
0.118 In. 


21 


2.71 In. 2.71 In. 2.71 2.71 In. 271 in. 2.71 In. 


A, = 1.9608 Sq In. 
T,=0.9204 In.* 

= 1.382 In. 
= 0.731 In. 


(a) SECTION 1 


A, = 1.9608 Sq In. 
1,=0.4914 in.* 
Zy2= 1.056 In. 


= 1.056 In. 
(b) SECTION 2 


A3= 1.6902 Sq In. 
I= 0.4895 In.* 

Z y3= 1.136 In. 

= 0.544 In. 


(c) SECTION 3 


Fig. 19.— VARIATIONS OF THE Ris SECTION 


the horizontal thrust Hx and the end moment Mx caused by a temperature 
rise At® can be determined. 

It is usual to consider the effect of shrinkage of the concrete as equivalent 
to a fall in temperature of a specified At®. Hence, the stress caused by shrink- 
age can also be determined by the same procedure. 

Differential Temperature Change Between Ribs and Shell.—If there is a tem- 
perature difference between the outside and inside air of a shell roof structure, 
the ribs, exposed mostly to the outside air, will have a different average tem- 
perature than the shell. The difference in temperature between the shell and 
the rib can be considered to be At®. Obviously, stresses will be produced in 
the structure. 

An analysis of this case resulted in thermal stresses as great as 100 lb per 
sq in. for a temperature difference of 10° F in a hangar of the type built at 
Rapid City (340-ft span). 

The actual temperature distribution will be somewhat different from the 
assumed one because the temperature will vary continuously. Temperature 
measurements on an actual structure should be made to determine the variation. 
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Stability of the Structure-——The thickness ¢t of the shell is mainly governed 
by stability considerations. For a certain spacing of the ribs, /, a radius r of 
the shell, and a given distributed load, a minimum thickness ¢ of the shell is 
required (13) (16). The resulting membrane stresses reach only a fraction of 
the allowable concrete stresses in compression. It is of interest to note that 
the ratio of the thickness to the radius of modern shell roofs is smaller than the 
corresponding ratio for an egg shell. This fact illustrates the perfection to 
which these structures have been developed. 

The ribs, reinforcing the shell at regular intervals 1, increase the buckling 
stiffness of the shell panel. Furthermore, the ribs are indispensable for carry- 
ing a one-sided live load (such as snow) acting on the structure. To provide 
sufficient buckling safety for the entire structure, the ribs must have a specific 
minimum stiffness. The buckling load is proportional to the moment of inertia 
of the effective section. Previous considerations concerning the stresses caused 
by a volume change led to a stiffness of the rib that was as small as possible. 
Obviously the stiffness should be kept at the minimum required for stability 
and bending strength so that the stresses caused by a volume change do not 
become excessive. 

Consequently, the factor of safety against buckling of the ribs is important. 
A smaller factor of safety will require less stiffness, and therefore smaller 
stresses will be set up by volume changes. Plastic flow of the concrete affects 
the geometrical shape of the structure, and thus the forces, in a dangerous 
manner if the factor of safety is decreased excessively. In addition, the sec- 
ondary moments caused by deformations of the structure can be correlated 
to the factor of safety against buckling (17) (18). 

Location of the Shell with Respect to the Rib.—The idea of placing the shell 
in the middle of the rib was advocated by C. 8. Whitney, M. ASCE, in 1950 
(19). Considering the stresses caused by a volume change, Mr. Whitney 
concluded that the stiffness of the ribs should be a minimum. By placing the 
shell in the middle of the rib the stiffness of the effective section is diminished 
to a great extent, as shown in Fig. 19. The same decrease can be obtained, 
however, by taking a smaller rib, with the shell at one edge of the rib. To 
illustrate this decrease, the sections shown in Fig. 19 must be considered. 
Sections 1 and 2 are identical except for the location of the rib. Placing the 
shell in the middle of the rib reduces the moment of inertia by 47%. Section 3 
has a moment of inertia approximately equal to section 2 as a result of the 
reduction in the height of the rib. The normal force and bending moment 
caused by the volume changes are proportional to the moment of inertia of the 
effective section. (Compare Eqs. 37, 39, and 43.) The fiber stresses, how- 
ever, are approximately proportional to the distance z of the fiber from the 
neutral axis (Eq. 44). The ratio of the corresponding stresses of the two sec- 
tions is thus equal to the ratio of the corresponding distances z. Therefore, 
for the upper fibers, 

Gus _ _ 1.136 


for the lower fiber, 
ors _ 213 _ 0.544 _ 
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Section 3 has in the upper fiber a stress 8% higher and in the lower fiber a 
stress 48% lower than the stresses in section 2. The objection that in section 3 
important cross bending stresses o, are set up can be rejected, for in the most 
extreme case go; is 1.73 times the lower fiber stress of the rib (146). Therefore, 


from which it can be seen that the cross bending stress ¢, is 10% smaller than 
o.». Furthermore, the cross bending is concentrated over a very short length. 

Therefore, for section 3 with a 21% lighter rib than section 2, the bending 
moment and normal force caused by volume changes are identical. Neglecting 
the influence of the normal force, the stresses of section 3 are +8% and —48% 
of the stresses of section 2. Cross bending stresses a, are 10% lower than the 
maximum stress in section 2 and of local importance. (0, has the form of a 
strongly damped oscillation starting from the rib.) Thus, from theoretical 
considerations, section 2 does not have any advantage over section 3. Con- 
trarily, 210% of the dead weight of the rib can be saved by using section 3. 
Therefore, it ean be concluded that it is not only ineffective but wasteful to 
place the shell at the center of the rib. 

CONCLUSIONS 

Theoretical solutions for some problems encountered in designing long-span 
cylindrical shell roofs were presented. Tests on a steel model in the scale of 
1 to 30 of an actual structure built of reinforced concrete (the hangar in Rapid 
City) confirmed the theoretical analysis. Possible objection to the fact that 
the test results were obtained from a perfectly elastic model (only elastic strains 
were induced in the model), whereas the material (reinforced concrete) in the 
actual structure exhibits quite different properties, can be answered as follows: 


It is common practice to compute the forces and moments in statically 
indeterminate concrete structures by assuming the concrete to be a perfectly 
elastic material. The moment of inertia for a cross section is taken for the 
uncracked concrete section. On the basis of the computed moments and 
normal forces, the stresses in the concrete and the reinforcing steel are deter- 
mined under the assumption that the conerete is ineffective in tension. This 
can be done by using either the ‘‘n-theory” or the “‘limit-design theory.” 

Therefore, the similitude between the theoretical direct forces and the 
bending moment in the steel model and the concrete structure can be made. 
These values were checked by results of the tests performed on the model, 
built of a nearly perfectly elastic material. Between the stresses in the model 
and in the actual structure no such direct relationships can be derived. It is 
an established fact (proved by many tests), however, that reinforced concrete 
structures, analyzed as elastic structures and reinforced according to the com- 
puted moments and normal forces, behave essentially as predicted. 

Many tests, not reported herein, were made during the course of the investi- 
gation. Similar procedures of analysis, based on the use of the effective width, 
were found to be adequate in checking test results to a close approximation in 
all cases. 

It was found satisfactory to estimate roughly the value of n, since the 
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effective width for the proportions of structures tested and the type of stress 
variation induced by usual loads put n in a region where variations caused little 
or no change in the effective width. Other engineering uncertainties are of 
much greater magnitude. In applying this analysis to different problems, it 
may be necessary to perform analyses by use of a Fourier series representation 
of the stress variation along the juncture of the rib and shell. In such cases 
effective widths and stresses could be determined for each of the significant 
terms of the Fourier series, and the actual stress at any point could be deter- 
mined by superposition. 
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APPENDIX I. NOTATION 


The following symbols adopted for use in the paper and for the guidance of 
discussers, conform essentially with “American Standard Symbols for Struc- 
tural Analysis’”’ (ASA-Z 10.8-1949), prepared by a committee of the American 
Standards Association with Society representation, and approved by the Asso- 
ciation in 1949: 


A = cross-sectional area; 
effective width ‘of the shell; 
modulus of elasticity ; 
total direct force defined in Eq. 5; 
= horizontal thrust at the center of a rib shown in Fig. 8(a); 
Hx = horizontal thrust on rib shown in Fig. 15(a); 
H'x = virtual horizontal thrust; 
moment of inertia; 
moment of inertia at the center of a rib; 
= a function of A and Bl, graphically represented in Fig. 4; 
span of the roof; 
overhang distance; 
= total bending moment: 
M’ = virtual bending moment; 
Ms, = bending moment on a unit; 
bending moment at the center of a rib shown in Fig. 8(a); 
restraining moment; 
= virtual restraining moment; 
bending moment of rib in statically determinate system; 
twisting moment; 
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total normal force: 
N’ = virtual normal force; 
N,. = normal force of rib in statically determinate system ; 
N, = direct force per unit width of shell in the axial direction ; 
Ng = direct force per unit width of shell in the circumferential 
direction ; 
Noz and Nzg = direct shear forces per unit width of shell in cir- 
cumferential and axial direction, respectively ; 
number of complete waves made by a harmonic function around the 
circumference of a cylinder; 


= virtual load; 


applied uniformly distributed load; 
radius of the cylinder; 
radius of the effective rib section; 
distance; 
shear force on a unit; 
eoordinate in axial direction of shell; 
distance of a fiber from the centroidal axis of the effective section; 
distance from the loads V4, to the centroid of the effective section; 
coefficient of thermal expansion ; 
shell constant defined in Eq. 4; 
known horizontal movement; 
known rotation of the abutment; 
6x1 = rotation of the left abutment; 
6x, = rotation of the right abutment; 


= coefficient of elastic restraint; 


parameter defined by Eq. 3; 

circumferential stress along the connecting line of rib and shell; 
o, = circumferential stress in the lower fiber; 
og, = cross bending stress; 


Tez and Tzg = shearing stresses; 


o= 
ox = 


(1) 
(2) 
(3) 


(4) 


angle defined in Figs. 5 and 8; 
angle defined in Figs. 5 and 8; and 
angle defined in Fig. 8(6). 
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